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The stability of multi-electron bubbles in liquid helium is investigated theoretically. We find that 
multi-electron bubbles are unstable against fission whenever the pressure is positive. It is shown 
that for moving bubbles the Bernoulli effect can result in a range of pressures over which the bubbles 
are stable. 
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I. INTRODUCTION 

Multi-electron bubbles in liquid helium were first ob- 
served by Volodin et al.^ In their experiment a layer of 
electrons was held in place just above the free surface 
of a bath of liquid helium by an electric field. The field 
was produced by a positive voltage applied to an elec- 
trode immersed in the liquid. The electrons remained 
outside the helium because for an electron to enter liquid 
helium it has to overcome a potential barrier of height 
approximately 1 eV.^ When the field reached a critical 
value, the surface of the liquid became unstable and a 
large number of electrons entered into the liquid through 
the formation of bubbles. Each of these bubbles typically 
contained 10^ ^ 10^ electrons. The multi-electron bub- 
bles are of interest because they could possibly provide a 
way to study a number of properties of an electron gas 
on a curved surface ^2, 

As a first approximation, one can consider that the 
radius of a spherical multi-electron bubble (MEB) is such 
as to minimize the sum of the energy associated with the 
Coulomb repulsion of the electrons and the surface energy 
of the bubble. This gives an equilibrium radius of 
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where Z is the number of electrons, a is the surface ten- 
sion of helium (0.36 erg cm~^ at 1.3 K)^ e is the di- 
electric constant (1.0573 at low temperature), and the 
applied pressure has for the moment been taken to be 
zero. Thus, for example, for N = 10^ the radius is 106 
^m. 

So far, there have been a very limited number of ex- 
perimental studies of these bubbles ii^i^ii In this paper 
we first consider the stability of an MEB that is at rest 
in the liquid (section II). We find that, at least when the 
simplest model of the energy of the electron system is 
used, the bubble is unstable against fission whenever the 
applied pressure is positive. In section III we investigate 
how the stability of a bubble is changed when it is moving 
through the liquid. We have been able to determine the 
region in the pressure-velocity plane where the bubble is 
stable. 



II. STABILITY OF BUBBLES AT REST 

Since MEB's were first observed, there have been 
several theoretical investigations of the stability of 
these objects. The first discussion was given by 
Shikin^ and further analysis has been given by Sa- 
lomaa and Williams,— and Tempere, Silvera and 
coworkers j "'-?''"3ii-^i In the simplest model, the electrons 
are taken to be distributed over the inner surface of the 
bubble in a way such that the electric field is everywhere 
exactly normal to the surface. This ensures that the 
charge distribution is in equilibrium. The electrons are 
treated classically and so are localized at the surface in a 
layer of zero thickness (see below) . Thus the total energy 
of the bubble is taken to be 

E^Es + Ev+Ec- (2) 

Here Es — oS is the surface energy with S the surface 
area and a the surface tension, Ey = PV the volume 
energy with P the applied pressure and V the bubble 
volume, and Ec is the Coulomb energy given by 

Ec^j'-^^V. (3) 

Since the electrons can move freely around the surface, 
the field 8 inside the bubble must be zero and so the 
integral in Eq. ([3]) can be restricted to the region outside 
the bubble. If the bubble is spherical, the bubble radius 
that gives the minimum value of the energy is the solution 
of the equation 

^'-\lQ^ea + %^ePRo) ' 

For zero applied pressure, this gives the total energy of 
an MEB as 

E^^-{2^Z^e^ae^f'\ (5) 

Since the energy is proportional to Z^/^ the energy is 
always reduced if the bubble breaks into two. Hence, in 
the discussion of stability given here we are not consider- 
ing whether the energy of the bubble can be lowered if it 
breaks into pieces, but are trying to determine whether 
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there is an energy barrier that prevents the bubble from 
breaking. 

To consider whether the spherical shape is stable, write 
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where ryz,-™ = It is straightforward to show that to 
second order in the parameters rjim, the three contribu- 
tions to the energy can be written as 
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Hence the total energy is 

167r 2^ I 

E^UnaRl + ^PRl + Y^-ai ^ jryz^p, (10) 
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where the spring coefficients a; are given by 



ao = 2,aRl + 2PRI, 
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and 



ai^{l-2){l-l)aRl~2{l-l)PRl, (12) 



for ^ > 1. From this one can see that the bubble is 
stable against spherically symmetric perturbations pro- 
vided that iaRl + 2PRl > 0. This leads to the condition 
P > Pr where 
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For Z = 1, the spring coefficient ai is zero; this is to be 
expected since a perturbation of the form rjimYimid, 4>) 
corresponds to a simple translation of the bubble in some 
direction. For I — 2 the spring constant a2 is zero if 
the pressure is zero, and so this analysis of the effect of 
small perturbations to the initial spherical shape does 
not determine the stability of the bubble. The higher 



I spring constants are all positive at zero pressure but 
each becomes negative if the pressure is increased to a 
sufficiently positive value. It was noted by Tempere et 
alJ^ that if the pressure is negative (but not negative 
with respect to Pc), all of the spring constants will be 
positivei^ and so the bubble must be stable. 

The stability of the bubble at zero pressure is of espe- 
cial importance since in the experiments that have been 
performed so far there has been no applied pressure apart 
from the very small hydrostatic pressure due to the dis- 
tance the bubble is below the free surface. At zero pres- 
sure 02 is zero, and so we need to go beyond the low- 
est order in perturbation theory in order to investigate 
the stability of an MEB at zero pressure. One approach 
would be to calculate the terms in the energy that are 
fourth order in the rjim parameters. Instead we have per- 
formed numerical calculations of the total energy as a 
function of bubble shape. 

To do this, we describe the shape of the bubble using 
Eq. ^ but now do not restrict the parameters rjim to 
being small. When the bubble shape changes, the elec- 
trons will redistribute themselves over the surface so as 
to minimize the energy and to make the electric field in- 
side the bubble zero. For each choice of shape we use the 
finite element methocU^ to compute the surface charge 
distribution and the Coulomb energy. The simulation 
uses 1280 triangle patches. We start with a spherical 
shape and vary the parameters rjim to see if a state of 
lower energy can be reached without passing over a bar- 
rier. We have done this using a maximum value of 5 for 
I in Eq. This process was then repeated for a series 
of different pressures. We also performed similar calcu- 
lations with a maximum value of / of 15 but taking only 
TO = 0. Both procedures gave the same results for the 
stability. 

The result of this investigation is that for all positive 
pressures there is no barrier to fission, whereas for nega- 
tive pressures there is a barrier. This result holds for all 
values oi Z. To illustrate the path to fission, we describe 
results obtained for a simplified calculation in which only 
I = and I — 2 contributions are retained. Thus we write 
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Within this simplified model, fission occurs when 02 — oq 
and the bubble develops a hole along the z-axis, i.e., takes 
on a donut shape, when 02 — — ao/3. In Fig. [T](a), we 
show examples of contour plots of the energy in the oq- 
a2 plane. The pressure is —0.03 mbar and Z = 10®. 
There is a stable minimum with a2 equal to zero, i.e., the 
bubble is spherical. When the pressure is zero (Fig.[T](b)), 
there is still a point in the plane at which the energy of 
the bubble is stationary with respect to both oq and 02 
(at Oo = 23.8 fim and 02 = 0), but it is now possible 
to reach the fission line from this point without passing 
over any energy barrier. Note that along this path there 
is, of course, an increase in the value of 02 but also a 
substantial decrease in ao- Once the pressure becomes 
positive (see, for example, Fig.^c)), there is no point in 
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FIG. 1; Contour lines of constant energy for a multi-electron bubble containing lO'' electrons for three different pressures. The 
energy spacing between contour lines is 0.05 eV. The energy is shown as a function of the parameters ao and a2 as defined in 
Eq. (O. The electrons are distributed over the bubble surface so as to minimize their energy. 



the ao-a2 plane where the energy is stationary. 



the equilibrium spherical shape is 



These results can be compared with the earlier calcu- 
lations by Tempere et alJ^ who also investigated the sta- 
bility against fission. They used an ingenious method in 
which the bubble was described by 6 parameters chosen 
so that the shape of a bubble undergoing fission could 
consist of two spheroids connected by a hyperboloidal 
neck. The choice of parameters was such that the bubble 
could vary from consisting of a single sphere, to an el- 
lipsoid, and then all the way to separated spheres. They 
minimized the total energy of the bubble by adjusting 
these parameters subject only to the constraint that the 
total length L of the bubble had to have a given value. 
They then investigated how the total energy varied with 
L starting from a value of L equal to 2Rq. If the en- 
ergy decreased monotonically as L increased from 2Rq 
to a large value, this indicated that the MEB was un- 
stable against fission. If the energy first increases before 
decreasing, this indicates that the bubble is stable. To 
simplify the calculation, Tempere et al. made the ap- 
proximation that the charge density was uniform over 
the surface of the bubble. They concluded that at zero 
pressure even though there is a mode of deformation (the 
I — 2 mode) which can grow without increasing the en- 
ergy of the bubble, there should be an energy barrier 
which prevents fission^ii whereas we find no barrier. 

This difference in the results arises from the treatment 
of the charge distribution on the bubble. If the bubble 
is assumed to have surface charge density that remains 
uniform when the shape changes, it is straightforward to 
show that the Coulomb energy for small changes from 
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In this case the spring constant a[ for the /-th mode (con- 
sidering only I > 1) becomes 



+ 2)aRl + 2PRl - 
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Comparing this with the spring constant a/ when the 
charge redistributes (Eq. (fT2|) ) gives 
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Thus for all modes, except I — and ^ = 1, making 
the approximation of a uniform surface charge gives an 
increase in stiffness and makes it harder for the bubble to 
undergo fission. The increase in stiffness is to be expected 
since a redistribution of surface charge can only lower the 
total energy. In Fig. [21 we show energy contour lines in 
the ao-a2 plane for an MEB with 10^ electrons at zero 
pressure calculated by taking a uniform surface charge. 
One can see that within this approximation the spherical 
bubble is stable. 

There are several physical effects that are not included 
in the simplified model used so far. It is possible that 
allowance for these effects would change the stability of 
an MEB at zero pressure. A more detailed consideration 
of the Coulomb energy (the total electron energy, to be 
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FIG. 2: Contour lines of constant energy for a multi-electron 
bubble containing 10^ electrons at zero pressure. The en- 
ergy spacing between contour lines is 0.05 eV. The energy is 
shown as a function of the parameters ao and 02 as defined 
in Eq. (|14p . The electrons are uniformly distributed over the 
surface area of the bubble. 



more precise) for a spherical bubble has been given by 
Salomaa and Wilhams^'^*' using the density functional 
formalism of Hohenberg and Kohn.^* This makes pos- 
sible the inclusion of the kinetic, exchange and correla- 
tion energies, but how these extra contributions affect the 
spring constants is not clear and is difficult to calculate. 
Salomaa and Williams show that these extra contribu- 
tions to the energy make a very small contribution to 
the energy when Z is large. For example, for Z = 10^ 
the extra terms make a contribution that is roughly 4000 
times smaller than the form Z^e^ /2eRo for the energy 
used in the simple model. The calculation could also 
be improved, for example, by using a density functional 
theory to treat the surface of the liquid helium, and by 
allowing for the penetration of the electron wave func- 
tion into the hquid. All of these effects appear to be very 
small corrections to the total energy and hence are un- 
likely to change the spring constants by a large amount. 
However, it is important to note that even a small cor- 
rection could lead to a positive value for a2 which would 
in turn lead to a finite (but small) energy barrier against 
fission. As an example, consider corrections that arise as 
a result of using a density-functional scheme to describe 
the helium. For a bubble with radius large compared to 
the thickness of the liquid-vapor interface the first cor- 
rection to the energy can be represented by considering 
the surface tension a to contain a correction that is pro- 
portional to the total curvature k of the surface. Based 
on a simple density functional scheme used previously, 
it is straightforward to show that the correction to the 



surface tension is Act = a'n, where 

a' = 0.9 X 10"^ erg cm'^ 
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and the sign of the correction is such that the surface 
tension is increased for a concave surface of the liquid. 
Inclusion of this term changes the total energy by an 
amount AE which for a bubble at zero pressure is given 
by 



AE 
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It is straightforward to show that the spring constant for 
an ^ = 2 deformation at zero pressure now becomes 



a2 = Ua'Ro. 



(20) 



Because a2 is now positive at zero pressure there will 
be a barrier against fission, but clearly for large Z (e.g. 
Z ^ 10*), this barrier will be very small. 



III. STABILITY OF MOVING BUBBLES 

The above results indicate that one way to stabilize 
an MEB is to produce it in liquid that is under a small 
negative pressure. We now consider an alternate way to 
maintain a stable bubble. A bubble moving through a 
liquid will be affected by the local pressure change asso- 
ciated with the liquid moving around it. For a spherical 
bubble moving at velocity 1; through an incompressible 
inviscid fluid with density p, the Bernoulli effect results in 
a pressure variation over the surface of the bubble which 
is given by^ 

P{e) =Po + lpv^{9cos^0-5) 
8 



(21) 



For a bubble in liquid that is at zero pressure far re- 
moved from the bubble {Pq = 0), this changes the shape 
of the bubble in two ways. The term proportional to 
Yoo by itself would provide a negative pressure around 
the surface of the bubble and since bubbles are stable at 
negative pressure, this contribution serves to stabilize the 
bubble. The second term gives a positive pressure at the 
poles of the bubble and a negative pressure around the 
waist. This pressure distribution will distort a spherical 
bubble so as to make the parameter 7720 in Eq. ([5]), or 02 
in Eq. (fT4)) . to be negative. This tends to stabilize the 
bubble since, as can be seen from Fig. [U for fission to 
occur a2 has to become positive. 

We have performed computer simulations in order to 
find the shape of moving bubbles and the range of ve- 
locity and pressure for which they are stable. We start 
with a guess at the bubble shape and then calculate the 
charge distribution on the surface. This then gives the 
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pressure APc\{0) exerted on the surface by the electrons. 
We then find the flow in the Hquid. To do this we expand 
the velocity potential as 



$(0) =^B,Pi(cos( 
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where Bi are some coefficients and the sum includes 
terms from Z = 1 to / = 20. The coefficients are deter- 
mined so as to give a velocity distribution in the liquid 
such that in the frame of reference of the moving bubble, 
the liquid velocity at the bubble surface in the direction 
normal to bubble surface is as close to zero as possible. 
This gives a pressure at the bubble surface of 



(23) 



where Pq is the pressure in the bulk liquid far removed 
from the bubble and APb(^) is the Bernoulli pressure. 
The net inward force acting on unit area of the bubble 
surface is then 



Po + APb{9) + 2(7 K - AP,i{9), 



(24) 



where k is the total curvature of the surface and APc\{d) 
is the outward pressure exerted by the electrons. Each 
part of the bubble surface is then moved inward a dis- 
tance proportional to this force, and the process repeated 
until the equilibrium shape is found. The calculation 
used a maximum value of I of 15. The calculation as 
just described is based on the assumption that the bub- 
ble shape and velocity field have axial symmetry around 
the direction in which the bubble is moving. In order to 
test this assumption, we also performed a calculation in 
which axial symmetry was not assumed. This calculation 
used values of Z up to 5 and |m| < 5. This calculation 
showed that the axially symmetric shape was stable. For 
an MEB with Z = 10® shapes for three velocities are 
shown in Fig. [3l We are able to perform the numeri- 
cal calculation until the bubble becomes concave at the 
poles. This is shown by the dashed line in Fig. [51 

In Fig. m the distance i?poio from the bubble center to 
the pole and the radius i?waist of the waist are shown as 
a function of the velocity. In Fig. [5] we show a plot of the 
region in the pressure-velocity plane in which the bubble 
is stable. This region is bounded by two lines. For small 
velocities there is a critical positive pressure at which 
the bubble undergoes fission. At negative pressures the 
bubble becomes unstable against expansion. For zero 
velocity this expansion is isotropic. 

Note that the change in the shape of the bubble even 
for a small velocity is surprisingly large. This comes 
about simply because the Bernoulli pressure contains a 
finite term varying with angle as Y2q{9, (f>) but the spring 
constant a2 for this pressure component is zero. Thus, 
for an MEB the changes in i?poie and i?waist are linearly 
proportional to the bubble velocity whereas for a gas bub- 
ble in a liquid the spring constant a2 is finite and so the 
changes in dimensions are proportional to the square of 
the velocity. 
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FIG. 3: The shape of a multi-electron bubble containing 10® 
electrons for bubble velocities of 3, 10 and 20 cm s~^. 
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FIG. 4: Distance i?poie from the bubble center to the poles 
and the radius -Rwaiat of the waist as a function of the bubble 
velocity. These results are for a bubble containing l(f elec- 
trons moving through hquid in which the pressure at large 
distance from the bubble is zero. 



The region of stability of bubbles containing a different 
number of electrons can be found by scaling the results 
shown in Fig. [5l The instability pressure PmstahiZ^v) 
can be written in the form 



(25) 



6 



0,10 



0,05 



^ 0,00 

H 

OS 

W -0,05 

CO 

W 
ed 

-O10 



-015 



FISSION REGION 



STABLE REGION 



_1_ 



_J_ 



EXPLOSION REGION 
J I I I 



10 15 20 25 

VELOCITY (cm s ') 



30 



35 



FIG. 5: Plot of the region in the pressure- velocity plane in 
which a MEB containing 10^ electrons is stable. The region 
is bounded by the lines on which the two different types of in- 
stability occur as described in the text. Along the dashed line 
the bubble becomes concave at the poles and the numerical 
calculations become inaccurate. 



where / is a dimensionlcss function, A = (cr'^e/e^)^/^ and 
B = p(eVcr'^e)i/3. Thus, for zero velocity the critical 
negative pressure at which a bubble becomes unstable is 
proportional to and at zero applied pressure the 

critical velocity at which the bubble becomes concave at 
the poles occurs is proportional to Z~^/^. 

We note that in this paper we have treated the liquid 
as inviscid although, of course, helium above the lambda 
point has a finite viscosity and below the lambda point 
the liquid still has a normal fluid component. At suffi- 
ciently low temperatures the density of the normal fluid 
becomes very small and, in addition, the mean free path 
of the excitations making up the normal fluid becomes 
comparable to the radius of an MEB. Under these condi- 
tions, it appears that the only effect of the normal fluid 
is to determine the mobility of an MEB and there should 
be no effect on the shape change or the stability. For a 
bubble with Z = 10^ the mean free path of the thermal 
excitations becomes equal to the radius at around 0.6 K 
and at this temperature the normal fluid density is less 



than the total density by a factor of 4 x 10~^. But as 
far as we are aware, there have been no experiments with 
MEB's at low temperature. 

At high temperatures where the helium is in the nor- 
mal state, the situation is not so clear. It is known that 
when the Reynolds number is large (but not so large that 
the flow becomes turbulent) the viscosity results in a thin 
boundary layer on the surface of the bubble and the pres- 
sure at the bubble surface is close to the value that would 
result from potential flow.^^ This general idea would sug- 
gest that the inviscid approximation should give reliable 
results for the stability of MEB's over a wide range of 
Reynolds number. To determine this range one could cal- 
culate the effect of viscosity using the method developed 
by Li and Yan^^ and applied by them to calculate the 
shape and drag on gas bubbles moving through a liquid. 
We have not attempted to do this. We note that, Al- 
brecht and Leiderer=^. in their experiments at 3.5 K found 
that the mobility of the MEB's was between one and two 
orders of magnitude larger than expected on the basis of 
ordinary hydrodynamics. The reason for this is unknown. 

In the experiments of Volodin et alX and Khaikinj^ 
which were performed at 1.3 K, velocities of the order of 
10'* cm s~* were reported. At these velocities the normal 
fluid component would be in the turbulent regime and 
the bubble is moving so fast that it should lose energy 
through the production of quantized vortex rings. 



IV. CONCLUSION 

We have examined the stability of multi-electron bub- 
bles in liquid helium and found that stationary bubbles 
at positive pressures are unstable. We show that because 
of the Bernoulli effect moving bubbles can be stable even 
at small positive pressures. 
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